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Pointlike objects cause many of the divergences that afflict physical theories. For instance, the 
gravitational binding energy of a point particle in Newtonian mechanics is infinite. In general 
relativity, the analog of a point particle is a black hole and the notion of binding energy must be 
replaced by quasilocal energy. The quasilocal energy (QLE) derived by York, and elaborated by 
Brown and York, is finite outside the horizon but it was not considered how to evaluate it inside the 
horizon. We present a prescription for finding the QLE inside a horizon, and show that it is finite at 
the singularity for a variety of types of black hole. The energy is typically concentrated just inside 
the horizon, not at the central singularity. 



I. INTRODUCTION 

It is a fundamental fact of general relativity (GR) that 
there is no such concept as the local energy of the gravita- 
tional field. The local effects of gravity can be removed by 
transforming to a freely falling frame. A neutral object 
at the origin of a freely falling frame will not experience 
any gravitational acceleration. Tidal forces will remain, 
but they only act on particles that are separated by some 
distance. Since gravity has no local effect, there exists no 
local energy. 

Various definitions of local energy densities can be 
made by making reference to special coordinate systems 
or background metrics. Heuristically, if we have an ob- 
server that we consider static we could use its acceleration 
as a measure of the local gravity. More sophisticated ap- 
proaches yield a variety of quantities that are useful for 
certain applications. The structure of GR is such that lo- 
cal quantities representing an energy do not exist, there- 
fore any attempt to define them must use use concepts 
that are not a natural part of the theory, i.e., a special 
coordinate system. Studying the asymptotic behavior of 
the metric, as done in post-Newtonian approximations, 
or the behavior of the metric at spatial or null infinity 
H II II leads to more useful and natural formulas for 
the energy. From these ideas, we are led to the idea of 
finding the energy inside a given finite region rather than 
the energy at a point. 

Quasilocal energy (QLE) is the energy inside a two- 
dimensional surface. The surface could be a sphere en- 
closing a star or a black hole, a small box enclosing some 
matter undergoing cosmological expansion, or a compli- 
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cated, even disconnected, surface in the spacetime. In 
this paper, we follow the method of Brown and York 
P, which derives an energy from a Hamilton- Jacobi 
argument involving the canonical action. This QLE has 
many useful properties. For example, it agrees with the 
Newtonian limit for a spherical star, is applicable in ther- 
modynamic problems |j|, Q , and the asymptotic limit at 
Euclidean infinity is the ADM expression for energy. Fur- 
thermore, it can be directly obtained from the Hamilto- 
nian of the same action principle (footnote 14 of 0) with- 
out the need for any other geometric structures. There 
are many formulas for other quasilocal energies |5( (and 
many references given therein), derived using different 
methods and often having different properties. 

We can define the quasilocal energy of the electric field 
in classical physics for comparison. The electric field of 
a point charge falls off with 1/r 2 , and the energy density 
equals the field strength squared. We have for the energy 
inside a surface of radius R 
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which makes an infinite contribution because the charge 
is pointlike. This problem remains in quantum electro- 
dynamics and requires renormalization, where another 
infinite quantity is subtracted to leave a finite remain- 
der. 

In GR, the situation is somewhat different because ob- 
jects of a given mass cannot have an arbitrarily small 
size. Once they become too small, they collapse to form 
a black hole and an event horizon forms. The resulting 
object is effectively the size of the event horizon, and 
outside observers are shielded from the infinities at the 
center. The QLE for a Schwarzschild black hole of mass 
M has the large distance limit 
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The QLE becomes undefined at R < 2M, the radius 
of the event horizon. The large distance limit suggests 
that the energy will diverge at the center, although we are 
protected from seeing this behavior by the event horizon. 
However, this is still something of a problem because an 
observer can fall in through the horizon in finite proper 
time and survive to see the interior of the black hole. It 
may be useful to have a definition of energy for observers 
inside the event horizon. 

A major issue that arises when defining QLE is which 
observers to use. Different observers, for instance ones 
that are accelerated rather than free- falling, will mea- 
sure a boosted energy. In this paper, we choose the ob- 
servers that are stationary with respect to the boundary, 
i.e., their four- velocity is perpendicular to the normal. 
These are the observers that will not stretch or squash 
the boundary, which would affect the energy. We will 
show that this prescription does not depend on the time 
slicing, and that the QLE so defined has useful properties. 
By comparing the QLE with the square root of the area of 
the boundary (which we refer to as the radius, although 
they are not equivalent), we can determine whether the 
radial coordinate is timelike or spacelike. This is simply 
the idea that if there is too much energy inside the cor- 
responding Schwarzschild radius, a horizon must form. 

The derivation in involved the boundary term of 
the Hilbert action. We write this action for a general 
spherically-symmetric and static metric and show how 
the QLE formula can be modified so that the bound- 
ary term is treated correctly in cither case. Surprisingly, 
when not coupled to other fields, the energy of the sin- 
gularity at the center of a black hole is zero. The energy 
climbs toward a maximum value at a radius inside the 
horizon, and at the horizon has an infinite downward 
slope. The charged black hole has a finite but negative 
energy at the singularity. 



II. THE BROWN- YORK QUASILOCAL 
ENERGY 

We now review the quasilocal energy defined in 0. 
The basic idea is the Hamilton- Jacobi method in classical 
mechanics of expressing the energy as a variation of the 
action with respect to the endpoints. The generalization 
to curved spacetimes results in the following definition 



E= - / d 2 x^/o z {k-k ) 
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where is the induced metric on the boundary 
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and a is its determinant. In the last equation is a 
future pointing timelike unit normal for the spacelike hy- 
persurface S whereas n M is an outward pointing spacelike 
normal to the boundary 3 B which is also normal to B if 



u ■ n — which is assumed in this definition. The con- 
stant k = 8irG is just a constant of proportionality, and 
in natural units is just 87r. 

The fc in the above equation is the trace of the ex- 
trinsic curvature of the two-boundary's embedding into 
the spacelike hypersurface S. The fco term is the energy 
of the vacuum, which must be subtracted to obtain the 
physical energy. In our case, we take flat space as the 
vacuum so that fco is the trace of the extrinsic curvature 
for the same two-boundary embedded in flat space. It 
is natural to add such a term because otherwise the in- 
trinsic geometry of the surface would contribute to the 
energy, even with no gravitational energypresent. 

E(R) has been computed already in pj for 4 dimen- 
sional spherically symmetric objects outside the event 
horizon. The metric of the spacelike slice can be writ- 
ten as 

ds 2 = f{r)- 2 dr 2 + r 2 (dn 2 ) (5) 
and the unit normal to constant r surfaces is 

n" = (0,/(r),0,0). (6) 
The extrinsic curvature of the two-boundary is 



-er"V Q n„ ; k 



(7) 



where the above covariant derivative is taken in the 
spacelike slice, and a is serving as a projection opera- 
tor; it can be found from 
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where is the metric of the spacelike slice 0. The 
only connection coefficients that we need are 



and we obtain 



(10) 



(ii) 



The reference term fco comes from setting f(r) = 1 so 
that we are embedding the sphere in flat space. For a 
spherical star with energy density p(r), we have |l6j | 



f(r) = y/1 - 2m(r)/r 

m(r) = 4-7T / r' 2 dr' p(r') 
Jo 



(12) 
(13) 



or for a black hole we simply have m(r) = M and f(r) 
1 - 2M/r. With 



r 2 sin ( 



the QLE becomes 

E{R) = R 
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as long as r is a spacelike coordinate. 

Specializing to the Schwarzschild case where m(r) — 
M, there is a horizon at r = 2M. Inside the horizon the 
r coordinate becomes timelike as revealed by inspecting 
There is a choice of which sign the unit normal 
should have inside the horizon; we choose 



= (0,-(2Af/r--l) 1 / 2 ,0,0) 



(16) 



which we will justify later by examining the boundary 
terms. The QLE becomes 



E(R) = R 



1/2" 
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inside the horizon. The QLE of the entire Schwarzschild 
metric is plotted in Fig. (JIJ and shows three striking fea- 
tures. 

First, the quasilocal energy at the singularity is zero. 
In Newtonian gravity, the energy of the gravitational field 
would diverge at the center for a point particle. So the 
nonlinearity of general relativity has removed this infin- 
ity, and gives us a picture where the singularity is not the 
most important feature of the black hole. We expected 
that the mass seen at infinity would reflect the mass of 
the singularity, but this seems not to be the case. The 
black hole looks like an extended object when we consider 
the second feature, that the QLE attains its maximum 
inside the horizon at a radius of 1 + ^-M. Most of the 
energy of the black hole seems to be "stored" just inside 
the horizon. 

The third striking fact is that the derivative of the 
QLE matches across the horizon, but is infinite there. 
The energy should be continuous on physical grounds, 
but the derivative might not be. In fact, if we had chosen 
the opposite sign of the normal, there would be a cusp 
in the QLE at the horizon. We will show that this is the 
correct choice in the next section. 

We note in passing that in three-dimensional spacetime 
the QLE is constant everywhere, so these features would 
be missing. In 2 + 1 spacetime, there is no horizon and 
the metric only possesses a conical defect at the center. 
This is because the Schwarzschild metric in n dimensions 



ds 2 



2m i , , 
1 ? dt 2 



2m 

y.71 — 3 



1 - — r ! dr ■ rdii 2 

In three dimensions with dil 2 = d9 2 and \fa = r the 
QLE is constant everywhere. 



III. RELATIONSHIP BETWEEN ACTION AND 

QLE 

The above calculation is not convincing because it is 
not clear that we have chosen the correct unit normal. A 
different choice would change the sign and give a much 
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FIG. 1: QLE computed inside and outside the event horizon 
for a Schwarzschild black hole. Both axes are in units of the 
mass M, and the horizon is at 2M. 



different result. To fix this ambiguity, we will go back 
to the derivation of QLE from the action. The extrinsic 
curvature formula for the QLE arises from the boundary 
term in the action. From the behavior of this term inside 
and outside the horizon, we can determine how to modify 
the QLE. 

We begin by considering a region of spacetime M . The 
spacetime is foliated by spacelike hypersurfaces E, which 
can play the role of moments of time. The region of 
spacetime will be bounded by hypersurfaces that we will 
simply refer to as t\ and £2 ■ Each £ has a spatial bound- 
ary B, which in this paper will always have the intrinsic 
geometry of a sphere. The time history of the boundaries 
B will be called the three-boundary 3 B. 

The three-boundary 3 B is the important boundary for 
the derivation of the quasilocal energy. In Brown and 
York's derivation, this boundary must be perpendicu- 
lar to the time slices E. Besides simplifying the cal- 
culation, this is also an important physical point: the 
observers whose velocities are normal to the time slice 
are the observers that will measure the QLE. These ob- 
servers should be at rest with respect to the boundary, 
and therefore the time slices should be perpendicular to 
the boundary. If they are not, then the observers will be 
boosted, and it is to be expected that they will measure 
a different value for the QLE. In this paper, we will relax 
the condition slightly but measure the same energy. 

To make clear how the boundary term gives rise to the 
QLE, we will restrict our attention to only a simple class 
of metrics. Spherical symmetry is imposed so that we 
can easily embed a sphere with a given surface area in 
the four-dimensional metric. The metrics we study are 
static, so that the only nonzero derivatives are radial. 
While we would like to generalize this derivation at a 
later date, this restricted version gives several interesting 
results and has the benefit of being easily understood. 

A general spherically symmetric and static metric can 
be put in the form 



ds 2 



-eN(r) 2 dt 2 + ef(r)- 2 dr 2 + r 2 (d9 2 + sm(0) 2 d(j) 2 ) 

(19) 
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where e is 1 outside the horizon and —1 inside. N and / 
will be chosen to be positive. 

To investigate the properties of the action, we will im- 
pose boundary conditions at some fixed r. The boundary 
term that we add is the one suitable to fix the induced 
metric on the three-boundary 3 B. The induced metric in 
this simple case is 

^jdx'dx 3 = -eN(r) 2 dt 2 + r 2 (d9 2 + sm{0) 2 d(f> 2 ). (20) 

The action, with a boundary term added to fix the 
metric on the boundary, is 



S = -L f d^x^—gR-- f d 3 xVW\® 

lK JM k JdM 

where the extrinsic curvature is 



and the trace simplifies to 



e = 7 y e y = - 



(21) 



(22) 



(23) 



The bulk term for the action using this ansatz for the 
metric is 



Sbulk = - I d l .r sin <-) 



(24) 



x [N/f - e(Nf + 2rNf + 2rfN' + r 2 f'N' + r 2 fN")] 
where primes denote r derivatives. The boundary term 



Sbt = - ( - J d 3 xsm6(-f(Nr 2 )'). 



(25) 



The boundary term can be converted into an integral 
over all 4 dimensions by also fixing the metric at r = 
and integrating the derivative of the boundary term. 

Sbt = f d 3 x(BT) = f d 3 xdr{BT)' + f d 3 x{BT) 

J J J r—0 

(26) 

Adding the two terms shows that in this case, the ac- 
tion has a very simple form. 



S= - I d 4 xsm(6>) 
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(j + ef)N + (2erf)N' 



(27) 



We have dropped the constant that comes from the sec- 
ond term in l|26|l . because it will not affect the final result. 
Varying N yields 



(2er/)'=(- + e/) 



(28) 



and substituting back into the action shows that the ac- 
tion is an integral of a total derivative. Also doing the 



angular integrations (trivial because of spherical symme- 
try) gives 



5=^ J(Ndt)(erf) 



(29) 



Following y| , we define the quasilocal energy as minus 
the second term in parentheses, so 

E(r) = -(erf) . (30) 

The values for the Schwarzschild metric are 

n(r) = f(r) = y/e{\ - 2M/r. (31) 

We also need to subtract the energy of flat space, which 
does not depend on e; we are embedding a sphere in flat 
space where there is no horizon. The subtraction term 
just has fco = — -. This will be the subtraction term 
used in the entire paper. The result for the QLE for any 
metric of the form considered is 

E(r) = r(l - e/(r)) (32) 

which reproduces the result in the previous section. 

IV. COORDINATE INDEPENDENCE 

We can relax the restriction on the form of the met- 
ric slightly, and consider what happens when the time 
coordinate is given an r dependence. The Brown- York 
derivation requires that the t and r coordinates be per- 
pendicular at the boundary, a condition which is violated 
by this transformation. Mo re g eneral derivations of the 
QLE have been considered |10L fill where this con- 
dition is eliminated. We will not consider this issue in 
depth, but simply use the transformation to show that 
our version of the QLE is not coordinate dependent, for 
coordinate transformations of this type. 

If we make the transformation t = t(T,r), then we 
can write dt — t^dT + t r dr, where subscripts denote 
derivatives. The metric becomes 

ds 2 = -eN 2 (t T dT + t r dr) 2 +ef- 2 dr + r 2 dfl 2 (33) 

= -e{t T N) 2 (dT+—dr) 2 +ef- 2 dr + r 2 dn 2 

written in a 3 + 1 form where the foliations are hyper- 
surfaces of constant r, which is appropriate for finding 
the induced metric on the three-boundary. Two good 
examples of this form are the Eddington-Finkelstein and 
Painlevc-Gullstrand coordinates. The action becomes 



S 



8tt 



(Nt T dT + Nt r dr){erf) 



(34) 



The integral is taken over a contour T which holds r fixed, 
so the final result is that the proper time Ndt has been 
transformed to NtxdT, the proper time written in the 



■5 



new coordinates. The QLE, —erf, is not changed. We 
will not relate this quantity to the extrinsic curvature 
because it is not necessary to our point here. 

The definition of the QLE is made of simple compo- 
nents. We embed a three-boundary in the space in such 
a way that at some moment of time, it has the intrinsic 
geometry of a sphere with a specified area. The geom- 
etry of the sphere does not change when moving along 
the time coordinate of the three-boundary. The QLE is 
the change in the action with proper time, both of which 
are invariants. The invariance of the result under this 
particular type of coordinate change is therefore not sur- 
prising. 



V. DESITTER SPACE AND BLACK HOLES 

The details of the derivation are not changed if a cos- 
mological constant is added. The term added to the ac- 
tion is 



Sec 



(35) 



which is proportional to N. When the equation of motion 
obtained by varying N is substituted back into the action, 
the same result as before is obtained but now f(r) is 
different. 

The metric is 



N(r) = f(r) = 



(36) 



There is a cosmological horizon at r = y a • ^ s ex P ec ted, 
the energy continually grows with increasing r. The hori- 
zon forms when the energy inside the surface grows larger 
than 



E(R) = -R. 



(37) 



This is the usual Schwarzschild radius expressed in a dif- 
ferent form. 

The deSitter-Schwarzschild solution has both a black 
hole horizon and a cosmological horizon. The black hole 
has a large amount of QLE inside a certain radius, but 
outside this radius the gravitational binding energy pro- 
vides enough shielding to bring the energy below c 4 R/G. 
The inner horizon forms at this radius. As one gets far- 
ther away, the cosmological constant begins to contribute 
noticeably. The cosmological horizon forms where the en- 
ergy has once again climbed above the necessary value. 



VI. REISSNER-NORDSTROM 

The Reissner-Nordstrom metric for a charged black 
hole is interesting because the metric behaves quite dif- 
ferently at the center from the uncharged case. The form 
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FIG. 2: QLE of deSitter space (positive cosmological con- 
stant). Both axes are in the same units, with a mass scale 
proportional to 1/ y/X. 
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FIG. 3: QLE of deSitter-Schwarzschild, a black hole in a 
spacetime with positive cosmological constant. The units are 
the same as the plot of the deSitter QLE, and the black hole 
has an unrealistically large mass so that details can be seen 
in the plot. 



of N and / is now 



N{r) = f(r) = ( 1 - — + C - 



2\ 1/2 



(38) 



where e is the charge of the black hole in natural units. 
There are now two horizons at r± = M ± y/M 2 — e 2 . 
The outer horizon is the same type as the Schwarzschild 
horizon. The inner horizon exchanges the signature of 
the t and r coordinates again, such that t is a timclike 
direction. The consequence is that the singularity is now 
avoidable. The inward radial direction is spacelike and 
not timclike, and so particles are not inexorably drawn 
into the singularity. 

Adding a new field into the theory will not change the 
definition of the QLE. The QLE only measures the grav- 
itational energy, and so only the gravitational action is 
important. Of course, the addition of a new field changes 
the metric. One of the most important characteristics of 
the gravitational energy is universality. All mass-energy 
contributes to gravity, and so the QLE measures the en- 
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FIG. 4: QLE of a Reissner-Nordstrom charged black hole. 
Both axes are in units of mass of the black hole, and the 
charge e 2 = 0.8M 2 

ergy of everything inside the surface (including the con- 
tributions purely from gravity). 

The striking feature of this case is that the energy be- 
comes negative within a certain radius. The QLE in ei- 
ther region where the time coordinate is timelike is 

+ (39) 

and so the energy becomes negative for R < e 2 /2M. 
This is always inside the inner horizon. The energy at 
the singularity is E(0) — — |e|. The singularity has the 
electric field of a point charge, and so using just classi- 
cal electromagnetism, the energy should diverge for small 
radius. However, the gravitational binding energy is neg- 
ative, and while the cancellation is not perfect it seems 
that the binding energy at least makes the energy at the 
center finite. 

The geodesies of massive neutral particles in the space- 
time offer a probe of the effects of negative gravitational 
energy. The radial geodesies obey the equation 

r 2 + V{r) =pl-l ; V(r) = (e 2 - 2Mr)/r 2 (40) 

where po is the conserved energy per unit mass of the par- 
ticle, and r is the proper time derivative of r. A particle 
that starts from r — oo with zero velocity will not reach 
the center, but turn around at r — e /2M. Particles with 
higher energies will penetrate farther toward the center, 
but massive particles of all energies are repelled. This 
result is well known |14| . The turnaround radius agrees 
with the radius where the quasilocal energy becomes neg- 
ative, so it seems that the two effects are very likely con- 
nected. Negative energy densities are expected to possess 
repulsive gravitational fields, and ne gativ e gravitational 
energy itself should be no exception. [lj| 



VII. CONCLUSION 

We have shown that there is a sensible way to extend 
the definition of the quasilocal energy to surfaces inside 
an event horizon. The Schwarzschild singularity has zero 
energy, and the energy of the black hole mostly resides in 
a region just inside the horizon. The addition of a posi- 
tive cosmological constant does not change these features 
but adds a cosmological horizon. In these cases, when 
the energy inside a given radius is less than r (in natural 
units), the space and time coordinates play their usual 
roles. When the energy exceeds this quantity, a horizon 
forms and space and time switch roles. The derivative of 
the energy with respect to r at a horizon always seems 
to be infinite. These two features make it easy to locate 
the horizons on a plot of the QLE. 

In a Reissner-Nordstrom black hole, the singularity 
at the center behaves like a point charge, and so there 
should be a divergence from the positive electric field en- 
ergy. However, the contribution from the gravitational 
binding energy is negative and apparently cancels the di- 
vergence, rendering a finite energy at the singularity. The 
QLE is negative inside the radius e /2M, which is al- 
ways inside the inner horizon of the Reissner-Nordstrom 
metric. A massive neutral particle released from rest at 
infinity will fall to the radius where the energy becomes 
negative, then reverse direction and be repelled. This 
provides the justification for the QLE's negativity, which 
is also related to the effect that clocks inside this radius 
run faster than those at asymptotic infinity. 

In this paper, we have used a specific preferred set of 
observers to define the quasilocal energy. A direction for 
future research is to remove this restriction as in other 
work [ToL IllL Il2t Il3| to define a more invariant quan- 
tity. We would like to extend this work to non-spherical 
boundaries and more general spacetimes. Two of interest 
are the spinning black hole metric and the metric for a 
star collapsing to form a black hole. There may also be 
applications to semiclassical and quantum gravity. Also, 
the issue of under what conditions the quasilocal energy 
is negative and what this means physically requires care- 
ful attention. 
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